Hexagonal vs. Rectilinear Grids for Explicit Finite Difference Schemes for the Two-dimensional Wave Equation

Finite difference schemes for the 2-D wave equation operating on hexagonal grids and the accompanying numerical dispersion properties have received little attention in comparison to schemes operating on rectilinear grids. This paper considers the hexagonal tiling of the wavenumber plane in order to show that the hexagonal grid is a more natural choice to emulate the isotropy of the Laplacian operator and the wave equation. Performance of the 7-point scheme on a hexagonal grid is better than previously reported as long as the correct stability limit and tiling of the wavenumber plane are taken into account. Numerical dispersion is analysed as a function of temporal frequency to demonstrate directional cutoff frequencies. A comparison to 9-point compact explicit schemes on rectilinear grids is presented using metrics relevant to acoustical simulation. It is shown that the 7-point hexagonal scheme has better computational efficiency than parameterised 9-point compact explicit rectilinear schemes. A novel multiply-free 7-point hexagonal scheme is introduced and the 4-point scheme on a honeycomb grid is discussed. [email protected]